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Abstract 

This paper describes four computer programs developed to assist 
logging engineers to plan transportation in a forest. The objective 
of these programs,, to be used together, is to find the shortest path. 
through a transportation network from a. point of departure to a 
destination. Three of the programs use the digitizing and plotting 
capabilities of a programable desk-top calculator system to conven- 
iently enter and store a numerical description of the network to be 
analyzed. The fourth program solves the problem, by an efficient 
optimization algorithm. The BASIC language programs are described 
through an example problem, and complete listings are provided. 

KEYWORDS: Road administration/planning (forest), transport system 
design (forest), computer programs/programing, computers 
(desk -top) . 
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Introduction 

The determination of the shortest 
route (or path) through a network of 
available routes is often an important 
step in planning transportation. 
Although problems in finding the 
shortest route frequently arise in 
other disciplines (Bradley 1975) , 
they are most commonly associated 
with road networks, and this paper 
will discuss these problems. 

Since the mid-1950 ? s, much 
attention has been devoted to problems 
of finding the shortest route. Recent 
publications have discussed procedures 
for solving such problems by hand 
(Mandt 1973, 1974) and by large, 
digital computers (Eisner et al. 
1975, Martin 1963, Schnelle 1972). 
The Forest Service has recently 
adopted a network analysis system 
(Sullivan 1974) , which is perhaps 
the most flexible system in existence, 
but it requires the availability of 
a very large computing facility. 

The programs discussed in this 
paper fall into the middle ground 
between procedures intended for hand 
computation and those developed for 
large-scale computers. These programs 
were developed and tested on a 
Hewlett-Packard 9830 programable desk- 
top calculator and were designed to 
take advantage of the digitizer- 
plotter capability of that calculating 
system. The shortest path program 
itself, however, is quite general and 
could be used, without significant 
revision, on any calculator or com- 
puter that uses the standard ASCII 
(American Standard Code for Information 
Interchange) BASIC programing language. 

Network Analysis 

A transportation system (such as 
a system of forest roads) may be 
described as a network, a collection 
of interconnected segments or links. 
Each link describes a unique path 
between two adjacent nodes. A node 



is any feature that might be treated 
as the point of departure or destination 
of some path through the network, such 
as a landing or mill. Nodes may also 
be road intersections, viewpoints, 
scaling stations, and bridges. Nodes 
are also commonly used to indicate 
points at which road design standards 
change or there is a marked change in 
grade or curvature. Such changes 
would be expected to influence costs 
of hauling logs (Byrne, et al. 1960) 
and may therefore be of interest in 
the solution of many forest transpor- 
tation problems. 

Figure 1 is a map of a small but 
representative road network for a 
hypothetical forest. Nodes have been 
superimposed on the network at each 
landing, at two grade breaks, at one 
point at which the road standard 
changes, at a bridge, at each of the 
three road intersections, and at the 
destination (the mill) . Numbering 
of these nodes is arbitrary. Each 
link, however, is referenced according 
to the adjacent nodes which it inter- 
connects. Furthermore, links are 
commonly bidirectional. That is, 
link 1-2 is "different' 1 from link 
2-1, even though each represents the 
same segment of road. This is, a useful 
distinction for differentiating the 
cost of hauling logs over an adverse 
grade in one direction and hauling them 
over a favorable grade in the opposite 
direction along the same link. 

SHORTEST-PATH PROBLEMS 

The specific problem addressed in 
this paper is that of finding the 
shortest path through a network from 
a specified point of departure. The 
unit of measure selected to judge a 
path f s length can be almost anything. 
Common units are hauling costs, distance, 
or time; but other units, such as con- 
struction cost, maintenance cost, or 
even a measure of the scenic or esthetic 
value along the link, could be used. 
Although this problem would normally 
be solved for the shortest path between 
a specified point of departure and a 
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Figure 1. Source map for the example problem. 



specified destination, it may occasion- 
ally be of interest to solve the prob- 
lem for all the shortest paths from the 
departure node to all other nodes. As 
an example, we might use this capability 
to solve for the best network of back- 
haul routes (i.e., when a truck is trav- 
eling empty) from the mill to all 
landings. The problem of finding all 
the shortest paths from a departure 
node to all other nodes is called the 
"n- shortest -paths problem/' It seems 
that this problem would be much harder 
to solve than the simple shortest path 
problem, and several of the faster 
algorithms for solving the shortest 
path problem also find the n shortest 
paths at the same time (Bradley 1975) . 
The solution procedure we have selected 
belongs to this class of algorithms. 
It is called the "Moore algorithm" 
and has been described fully by Martin 
(1963) . Although it always detennines 
the n shortest paths from a specified 
point of departure, our program per- 
mits the user to decide whether all 



Network Analysis Programs 

The programs (to be used together) 
are of two types : three programs are 
of the first type; they essentially set 
up the problem by generating the network. 
The second type (the fourth program, 
"Optimization Program") solves the 
problem. The network generator programs 
are machine dependent; that is, they 
were developed specifically for the 
digitizing, plotting, and computing 
capabilities of the Hewlett-Packard 
9830 system. Numerous other systems, 
however, have similar capabilities; 
the description in this paper is 
intended to be general enough that 
the programs can be converted to other 
systems . 

PROGRAM LIMITATIONS 

In their present form the programs 
can accommodate a maximum of 60 nodes 
and 255 links, with no more than 8 links 



limitations are a function of the 
calculator memory and could be increased 
if memory capacity were expanded. The 
calculator systems on which the pro- 
grams have been developed and tested 
have the following configuration: 

HP-9830 calculator with 15,808 
bytes (7,904 words) of read-write 
memory 

Additional read-only memories: 
plotter control 
extended input/ output 
string variable 

Flatbed plotter (9862A) 

Digitizer (9864A) 

Thermal page printer (9866A) 

NETWORK GENERATOR 

The most burdensome phase in 
network analysis by computer is 
usually data entry. Nodes and links 
must be identified and link values 
must be assigned before the analysis 
can begin. A primary advantage of 
the approach reported here is that 
the network definition task has been 
automated by the use of the HP-9830 
digitizer. 

Three network generator programs 
are provided for entering problem 
data into the calculator: one is 
loaded into the main memory of the 
calculator, and the other two are 
loaded into special function keys 
fp and , . Table 1 summarizes the 
input procedure for the data required 
when these programs are executed. 
Note that after all three programs 
have been loaded into the calculator, 
the program in main memory is 
executed first. This provides the 
user with some preliminary instruc- 
tions and initializes certain 
internal memory locations. The pro- 
gram on special function key fn, 
which enters the network via tne 
digitizer and reproduces it on the 
plotter, is executed next. This 
program also computes and applies a 
"rotation factor" so that the user 
need not be concerned about whether 



the source map is correctly alined on 
the digitizer. After the network has 
been completely digitized, the program 
on special function key f, is executed. 
This program establishes a "link matrix" 
which identifies individual links by 
their "from" and "to" nodes. It also 
prints a link summary and permits the 
user to store the link matrix on a 
cassette tape if desired. 

Figures 2, 3, and 4 summarize the 
printed output resulting from running 
the three programs listed above. This 
output is for the example program in 
figure 1, entered as shown in table 1. 
It is important to note that these 
tables contain printed output only. 
Numerous requests for data and certain 
instructions to the user are displayed 
(in the calculator display window) while 
these programs are running but are not 
printed. This convention produces a 
"clean" printout. A complete track 
of all printed and displayed messages, 
as well as a record of all inputs, can 
be obtained if desired by putting the 
calculator into "PRT ALL" mode before 
running the program. 

In addition to the printed output, 
the program on key fg also plots the 
digitized network on the flatbed 
plotter (fig. 5). This plot is useful 
as a visual check to insure that the 
network has been correctly entered into 
the calculator. 

SPECIAL CONSIDERATIONS 

It should be evident from the 
example (fig. 1) that it is often 
necessary to re-digitize nodes which 
have already been digitized in order 
to correctly describe all routes of 
the network. As an example, consider 
nodes 4 and 5. Node 5 must be digitized 
after node 4, and then node 4 must be 
redigitized in order to continue the 
network. When this happens, the 
calculator is programed to sense that 
the node has already been digitized, 
and a message to that effect is printed 
(fig. 3). 



Table \--lnput explanation for network generator programs 



Visual prompt on display or printer 
(oral prompt by digitizer) 


Keyboard or 
digitizer response 


Explanation 



DOCUMENT HORIZ.DIM. (INCHES)? 
DOCUMENT VERT. DIM. (INCHES)? 

SELECT FUNCTION KEY (FO) 



SET THE ORIGIN AT POINT 0,0 ON THE 
NETWORK SOURCE DOCUMENT 
(digitizer "beep") 

DIGITIZE POINT 8,5 
("beep") 



SCRATCH A. 
LOAD 5. 
LOADKEY 4. 
RUN. 



5. 



Depress red "0" 
button on cursor. 



Digitize upper 
right-hand corner 
of map. 



Digitize node 1. 
Digitize node 2. 
Digitize node 3. 
Digitize node 4. 
Digitize node 5. 
Redigitize node 4. 
Digitize node 6. 
Redigitize node 2. 
Redigitize node 6. 
Digitize node 7. 
Digitize map origin. 



Redigitize node 4. 
Digitize node 8. 
Digitize node 9. 
Digitize node 10. 
Digitize node 11. 
STOP. 



. . . (link matrix summary is printed) 
LINK MATRIX COMPLETE. STORE? YES. 



START DIGITIZING 


NODE LOCATIONS 


BEGINNING WITH 


YOUR CHOICE FOR 


NODE #1 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 




("beep") 





FILE # (MINIMUM SIZE =768)? 



0. 



Clear calculator memory. 
Load programs from tape. 

Begin program execution. 
Enter map dimensions. 



Depress special function 
key f to digitize network. 



Set point to which nodes will 
be referenced. 

Use "S" button on cursor. 



Use "S" button on cursor 
to enter node locations. 



Define link between node 
2 and node 6. 

Use "S" button, not '0'; this 
informs the calculator that 
a new route is to be entered. 

Enter new route, 
starting at node 4. 



Halt digitizer program. 
Depress special function key 

fi to print a link assignment 

summary. 



Indicate that the link matrix 
is to be stored on tape. 

Store link matrix on tape 
file 0. 



. . . (program terminates) 
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Figure 2. Printed output from the program in main memory (loaded from 
tape file 5) for the example problem. 



SET THE ORIGIN FIT POINT 0s. 6 ON THE NETWORK 
SOURCE DOCUMENT NITH THE DIGITIZER,, 



DIGITIZE POINT 3 ! 5 

STFlRT HI GIT I ZING NODE LOG FIT I ON 8 BEGINNING NITH YOUR CHOICE FOR HO HI 



NODE FIT 


1 07? 


3,, 70 


NODE FiT 


2 n 7' 8 !' 


2 '4 5 


NODE FT!" 


3 9 6 i> 


-'"' "!' - ".' 


NODE FIT 


'"" C! 1 '"!' 


2 7 9 


NODE FIT 


7 i' 5 i. 


1. 72 


NODE FiT 


;::. 43 


I ,,42 



NODE 
NODE 
NODE 



FIT 
FiT 
FiT 



'::;' @ !' 
6 3 !' 
, 54s 



IDENTIFIED FlS # 

IDENTIFIED FlS # 

I HEN TIE I EH FlS t 

I HE NT I FIE I) FlS # 

I HE NT IF I ED FlS # 



I HE NT IF I ED FlS # 
U,,4S il I HE NT IF I EH FlS tt 



4 ,,45 
4 6 1 
4 1 



RE COR. D EH EFlRLIER FlS # 



r.i i;:; .t.j. 

I I-:::' -It 
FlS tt 



RE COR HE I) EFlRLIER 
RE CO RUE I) EFlRLIER 

STFlRT NEW ROUTE, 

EH:- GIN RESTflRTS FiT 

THFiT WHS IDENTIFIED 



RECORDED EFlRLIER MS 



I DENT I FIE I) FlS # 8 

I HE NT IF I ED FlS tt 9 

I HE NT IF I ED FlS tt 10 

T T'I i:::' i...i "i" T ir T r: 1 'n i'"i C f :'' 1 ' 




:. % .:,:. # L J | 1 1 

FROM I IN I 



| E j-.j "|" '"; | | |v| |v| 

TO HI Mil: 1 1 



UE1 



i ::: i I.... U E C, 



I 1 ' 
'U 



9 
q 

1 
1 6 
1 1 



b 
4 



4 
9 

1 

9 
i i 
1 



(These columns are available as scratch 
space. They are intended for the user 
to write in values, such as cost, 
distance, or time, which may be assoc- 
iated with each link.) 



Figure 4. Printed output from the program on special function key f 
(loaded from tape file 4) for the example problem. 
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Occasionally it may be desirable 
to enter an entirely new route to 
avoid redigitizing several points. 
As an example, refer to the dead-end 
route leading from node 4 to node 11 
(fig. 1) . Rather than digitizing 
nodes 4, 8, 9, 10, and 11 and then 
retracing the route back to node 4, 
it is convenient to enter those nodes 
as a new route. This is done after 
all other nodes have been digitized. 
For exajnple, the last node to be 
entered is node 7 (figs. 1 and 3). 
Then, the fact that a new route is 
to be entered is signaled by redigit- 
izing the map origin (lower left-hand 
corner) , followed by the nodes of 
the new route (4, 8, 9, 10, and 11). 
Any number of new routes can be 
entered in this manner. 



OPTIMIZATION PROGRAM 

After the network has been fully 
described, the calculator memory is 
cleared and the optimization program 
is loaded into main memory to solve 
the problem. The link matrix created 
with the network generator may then 
be loaded from a tape (if it was 
stored as indicated at the bottom of 
table 1). Otherwise, it can be entered 
at the keyboard. Although the latter 
capability may appear unnecessary, it 
has been included to permit users 
without digitizers to use the optimiza- 
tion program. Such users would simply 
enter the link matrix at the keyboard 
without having previously run the 
network generator programs. 

Once the link matrix has been 
entered (either from a tape or through 
the keyboard) , values must be associated 
with each link. These values can also 
be entered from a tape or through the 
keyboard. Link values might have 
been stored on tape after a previous 
run of the optimization program; 
normally, they will be entered at 
the keyboard. As shown in table 2, 
the calculator rrroTnnts the user bv 



user is given the opportunity to cha 
any of the link values, to add links 
or to delete links. Either the link 
matrix or the value array or both ma 
then be stored on tape if desired 
(see table 2) . 

The optimization routine begins 
after the user identifies the initia 
node, or point of departure (table 2 
As discussed earlier, this routine 
finds the shortest paths from the 
initial node to all other nodes. Th 
method used for the optimization is 
not discussed in this paper; interes 
readers are referred to a complete 
discussion of the Moore algorithm 
(Martin 1963) . 

When the optimization is comple 
the calculator requests the number 
of the terminal node (see table 2) . 
The user may then indicate a specifi 
node, in which case the shortest pat 
from the initial node to that termin 
node will be priiited; or he may requ 
that the shortest paths to all nodes 
be printed. 

After these results have been 
summarized, the user is permitted to 
run a new optimization from a new 
initial node. He may also enter a 
new set of values, either from a tap 
or through the keyboard. 

Figure 6 is a listing of the 
printed output from the optimization 
program for the example problem, Fo 
the minimum path analysis, node 7 
(a landing) was selected as the poin 
of departure and node 1 (the mill) a 
the destination. The ''optimum path" 
from node 7 to node 1 is via nodes 
6, 4, 3, and 2 (see fig. 7). The 
total estimated tf cost n of following 
this path is 19. 



Table 2 Input explanation for the shortest path program 



Visual prompt on display or printer 


Keyboard 
response 


Explanation 



LOAD LINK MATRIX OFF TAPE? 

FILE # OF LINK MATRIX STORAGE? 
LOAD VALUES OFF TAPE? 



VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 
VALUE 



FOR LINK #1? 
FOR LINK #2? 
FOR LINK #3? 
FOR LINK #4? 
FOR LINK #5? 
FOR LINK #6? 
FOR LINK #7? 
FOR LINK #8? 
FOR LINK #9? 
FOR LINK #10? 
FOR LINK #11? 
FOR LINK #12? 
FOR LINK #13? 
FOR LINK #14? 
FOR LINK #15? 
FOR LINK #16? 
FOR LINK #17? 
FOR LINK #18? 
FOR LINK #19? 
FOR LINK #20? 
FOR LINK #21? 
FOR LINK #22? 



WANT TO CHANGE ANY LINK VALUES? 



WANT TO ADD A LINK? 
WANT TO ERASE A LINK? 



STORE LINK MATRIX? 
STORE VALUES? 

FILE # (NOTE: SIZE MUST BE >512)? 

*******MINIMUM PATH ANALYSIS 
INITIAL NODE IN PATH? 

TERMINAL NODE (999=ALL NODES)? 



SCRATCH A 
LOAD 6 
RUN 

YES 




NO 

3 
6 
9 
3 
4 
5 
3 
2 
8 
7 
2 
1 
6 
5 
5 
4 
3 
3 
2 
4 
1 
2 

NO 



NO 
NO 



NO 

YES 

1 



Clear calculator memory. 
Load program from tape. 
Begin program execution. 

The link matrix was previously stored 
on tape; prepare to retrieve it. 

Identify tape file number. 

The link values have not previously 
been stored. 

Enter values corresponding to the links 
as numbered in figure 4. 

Although the link values used in this 
example are whole numbers, decimal 
numbers may also be entered. 

Note that the values used here 
correspond to those in figure 1. 



This option is provided for correcting 
wrong entries; in this example, no 
corrections are necessary. 

These options are provided for making 
corrections or revising the network; 
here, no corrections or revisions are 
necessary. 

The link matrix has not been changed, 
and it is already stored on file 0. 

Indicate that the link values are to 
be stored. 

Store link values on file 1. 

Node 7 (a landing) is selected as the 
point of departure. 



1 Node 1 (the mill) is selected as the 

destination; if 999 had been entered, 
the minimum paths from node 7 to all 
nodes would have been printed. 

. . . (the minimum path summary is printed here). . . . 



LINK: 



JMMRRY FOR SYSTEM WITH 11 
FROM NO DEM TO HODEtt 



1 



NO I H 



! ! ! I I 



LINKS 



4 
4 



MINIMUM PRTH RHRLYSIS,. 



1 
ci 

1 1 
1 @ 



MINIMUM PRTH FROM NODE 7 TO NODE 
THE NODE TO NODE PFiTH IS 7 

6 

4 



19 



Figure 6. Printed output from the shortest path program (loaded from 
tape file 6) for the example problem. 
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Concluding Remarks 

The network analysis package 
described in this paper is intended 
for use with the small- to intermediate- 
size networks commonly encountered in 
planning transportation in forests. 
It is not suited to the analysis of 
large networks. Probably the most 
attractive feature of this package 
is the network generation capability 
which permits rapid, dependable entry 
of the network via a digitizer/plotter 
system. Also attractive is the 
interactive, problem- solving feature 
of the programs, made possible by the 
stand-alone capability of the desk-top 
programable calculator system. 
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Appendix 



Listing of the Network 
Generator Initialization Program 



TUN NC 2 "J? NIL S@!> 10 Is LIE 2 56s 3 ] 

REN NET WORK BUILDER. 

DIM fl*C30 3 

P R 1 N T T Fl 13 5 ' ' T H I S S E "I" F F :I R G R R M S 1 S H E S I G H E I) "1 ' ( j D 1 G I T I Z E Pi S E T F N D E ' ' 

:::i R I H "I" T ! ::: l EI 4 ' ' L C fl T I H S F F Fl S U R C E D C U i v l E N T i ::: ! S E i :; i C H I S D I G I T I Z E D I "I" I S ' ' 

F 1 R I H T T R 13 4 ' ' F 1 L T "1" E D ? N U II B E R E D i- Fl H D S T R E D I H T H E C i ::: l L C U L Fl T R . T H I S I S ' ' 

PRINT TRB4"fli::COI v IPLlSHED WITH Fl PROGRflM ON FUNCTION KEY Fi/i " 

P R I H T T ! ::: ! B'5 ' ' fl F T E R T H E N D E R E C R D I N G 1 S M P L E T E ? L I H K H U H B E R I H G I S D H E B 'i 
P R I H T T F! B 4 ' ' fl P R G R F! M S T R E. I) N F U N C "I" I N K E V ! :: 1 ' ' T fl B 1 @ @ 

P R 1 N T "I" R B 5 ' ' T H E F I R S T S T E P I S T I H F R M T H E F :I !.... T T E R F "I" H F : ::: : n U R C E D C U i v i E h 

@ PRINT T FIB 4 "SIZE " T Fl B 1 8 

MI: i ::N v ii::2 >e 

IUSP "DOCUMEHT HORZ I) I M (INCHES) 11 ; 

INPUT Gil 

DISP "DOC LINE NT VERT,, DIM,, (INCHES)"!! 

INPUT G2 

SCR LIE 8 n 1:1*10 On @i> G 2*100 

P R I N T T R B 5 ' ' W Fl I T F R T H E S E L E C T F U H C T I N K E Y ( ! :: ) N E S S R G E :> T H E N P R C E E D ' 

PRINT T FIB 100 

R E M * * * * * * ' ':''" * * I N I T I fl L I Z E T H IE N D E M Fl T R I X 

EOR 1=1 TO 60 

F' R J :;::: 3 T 1 

HI: i!.j:]~0 

NEXT J 

i0 NEXT I 

:0 DISP "SELECT FUNCTION KEY (PCD" !'. 

'0 STOP 

iO END 



Listing of the Network Generator Program 
on Special Function Key f 



,| i I 'i, 1 | 1 1 i i i i( i i ' ii i - ! ., i ii i 1 1 Miil! II! I 'I i 1 i I i i ! ' ' ' 

I'M i ; ! ;: i MiMV .iinl-'i \ I'm 'Nil Hi !! I 1 1: I III i ! i R ! ' ' - 
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', 'ii i -| M i i i I ill 1 - H "I 'I i, I 1 1 I- 1 1 ' i MRT3 FIT R NODE :! 

! !M i 'r'i I I i i I II ! i ' " i I ii i I Mil 1 ' i Ml:::. N i 1 i"' I !:::..! J i:::. Hh!l i. i:.;. 

1 , i ii i 1 1 ' i i i i i ii . i 



' | i I II i , I ! ' , i i | ! ! ! | i I 

'I i'l II I! ill! i' I I III M l.'ll 

'!,"i i '!' ill! i in I'l 1 i"i-'l i ill- 1 Ml H i l iRLIER FlS tt l! J:L 

I :u l ' I! 

Mil I- 1 1 i! ' ll I 1 1 1 M| PI 1 1 I M ,.' 1 1 1 1 , R Eh R IN ING F'O I UTS 
'I'M: i MI<' i Ml I I MI I linn 

-I'M Ml- 1 | | | i ' I ! ' > 

i ,'ii i 1 1 1 1 i' 1 - 'i- i ; ; i -I v i 



i F3,, @ 



: 1 @ 
S 2 @ 

) .::' '::!' 
549 

550 

: i f> 



530 
590 

5 

::- 1 



540 
350 

360 
578 
586 
i3 9 
700 
?' 1 
7 2 
i :: ' 3 13 
740 
7 5 
769 
770 
7 3 

7 9 
308 
8 1 
3 2 9 
338 

8 4 9 
3 5 8 
8 6 
8 7 
389 
398 
980 



R E l v i C H I:::. C l< T S E W H E T H E R R 1 G; I N H ! ::; i 8 

T F ft R R ' X 1 * X 1 !- Y 1 * Y 1 ) > -, 1 T hi E N 5 5 

PRINT TfiB45"STfiRT NEW ROUTE,, " 

GOTO 360 

Ji=M[ 1 3+1 

GO SUB 348 

IF J1=M[ 1 3+1 THEN 610 

i : -' p T H "I" T R E: 4 5 ' ' R E C R D E I) E i ::: i R !.... 1 E R i ::; l S ; i : i : ' : . J 1 

MC 1 3 = M[ 1 3 1 

GOT 3 670 

NRITE < 15 * 268>X1 ? Yl ;< Jl 

N[ Jl is 1 3~'X1#108 

NC Jl ! 2 >Y1#10@ 

PLOT Hi: Jl 9 1 IhHL J1-. 23 :0 

CPLOT 1 ii 1 

L.RBEL. OiO Jl 

PFI1 ESTflBLISH WHETHER THE I 1 HOD I MI|.' 

POP J=3 TO lO 

IF HC Ji ii J::KI:I THEN 790 

IF H [ J 1 I. J 3*1=0 THEN 748 

REN THE 11 NODE MRS NOT BEEN RE COR DEL 

HI: Ji I- J>n 

GOTO 778 

NEXT J 

P R I H T ' ' "I" N fl N Y !..., I N K 8 E N I N i ::: i T I H G F R i v l 

END 

PLOT HC II ii 1 In Hi: I In 23 ii 2 

PLOT N[ Jin 1 3i. H[ Jin 2] ii 1 

I1-J1 

Mi: 1 3=Mi: 1 3+1 

NEXT i< 

PRINT "NUMBER OF NO HE 8 EM I I H 1 . I ! I'll I i 




FOR Jl-1 TO M[ 1 3 

IF 8GRC (Hi: Jl ii 1 3 XI* 100) IV i > III I I , ' I 

PEN IF NODE WR8 NOT I DEN I I I' II H I. Ill- 1 

NEXT Jl 

RETURN 

END 
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Listing of the Link Matrix Generator Program 
on Special Function Key f- 
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i i ! Ill I I -. I HIM Mil II ,:-iM 

1 I ' "'' "'"I I' i '"'I i 1 ' 'i ' M|| Nil ,1" NOT HFlVE! BE:! EN BEEN RECORDED ,. NE DO S0 

i i ,i i ! i i , ' ,'i i 

i i' ' rl Mil II VI I -".I i Mil Mil 1 . REOORlJEIi PREVIOUSLY, 
i Ill ! I . i i ! 

iilll.il i i Mill i I i 
i : ! i I . ! I iil I l . I I 
i iii i 'i 'i iil I l -I , ! i 
i I ! M 1 1 I i I i !: 
i 1 i i ' I | i 

i MI ,; i M 

i Ii ! -I ' "I i iil Ni i i i- I " i i H'll i I IE:!,, STORE" S! 

! liM'lM III 

i ! i lii "VI >" i III II Mill 

I , i III ' 

i i 1 ! H I'M I'lll'l in M Ml- 1 ! 

i I i i'l Mil! 

i MI,!' "I 1 1 I ii i H M| I l-llll-l i 1 1 i SIZE^SER) " S 
i ! ili'NI I '' 

iii ' i i . i I 1 I 
ill 'i,i > l N! I I 

1 ' I i Mil, | IN I i 

i < i l I I ' ' i I I ! I i I I I ' ' 'i I 

i M NliiH 
i ! ilH 



Listing of the Shortest-path Program 



! HI i|' i I l-'i it ii''i iN. i i .1 h ii- 1 1 . i I- 1 1 M i i > ii i : M i 
! I' ."n... . I VM ,"!. I- Ml i.M hi M : ,| t.li |. I! | | ,,n ,, ,1 , , 
MI Pi NMMi 
4fi F 8 ::::: i r 9 ::::: 1 9 ;::;: L 9 ::::: 

58 DISP "LORI) LINK MR IRIX n| I M-ll i : "5 
66 INPUT Fl* 

70 IF R:ili-"NO" THEN 346 

88 DISP "FILE* OF LINK MFlTRIX STORflGE" 5 
90 INPUT F9 
1i"ii"i LOFlIi DlTlTl F9sL 
HO Ll^LC 256s 1 ] ' 
120 Nl :: Li: 2 5 Si' 2 ] 

130 DISP "LORD VRUJES OFF' TFlPE"! 
140 INPUT fit- 
150 IF R::N"HO" THEN 220 

1 i ::: :0 DISP "FILEtt OF VFILUE STORRGE" !i 
170 INPUT F8 

188 IF F3>F9 THEN 200 
190 RENIND 

2 L i ;:; l U D Fl T i ::; i F 8 ;> V 
2 1 I 9 ::::; i 

2 2 F 1 P. I H T ' ' * -i'i- ;; L I H l< S U N M Fl R V F R S V S "I" E II W I T H ' ' N 1 ' ' N ( ) D i:::. 3 ! ;:; 1 N D ' ' L 1 ' ' i.... 1 H I : 
230 PRINT "LINK* FROM NODE* TO NOD Lit ;; VFiLUE" 

2 4 F R l< :::;: 1 "I" L 1 

250 i i<[NTa,i: Kn 3 ::i/:i,oo) 

260 T1--M.J: Ki 3 ] 11*100 

270 IF 19^1 THEN 300 

280 DISP "v'Fll.JJE FOR LIHK#"K!! 

290 INPUT VCI<] 

300 PRINT K.H I 1 iiT'1 n V[ l< ] 

310 NEXT K 

320 IF L9tt0 THEN 590 

330 GOTO 5 50 

340 PRINT "KEVBOflR'D ENTRY NODE,, " 

350 DISP "NUMBER OF NODES" : i 

:~!60 INPUT Nl 

370 DISP "N UNBER OF LINKS"; 

330 INPUT LI 

: :i: 9 P R I H T ' ' * * # L I H K S U N 11 Fl R V F R 8 V S T E M N I T H ' ' i ! 1 ' ' N u J J 

1- 8 P R I N T ' ' L I N K tt F R l v l N D E 1* T H 1) E i : ! : s 

I 10 !....[ 25 6 :i 1 ]-l....l 

420 Li: 256;' 2 1 ::::: N1 

430 FOR IO1 TO LI 

440 DISP " INITIRL NODE FOR LINKir'Ki! 

450 INPUT II 
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, ,ii ! i MM 

,' in ni .i 1 ' i IN 1 1 111 mill 1 , i n|. 11 m n"i i ; i 

,' .M Mil Ml i I 

.'MI h i M " I I I 'I'll Mill III. !i| Mi! 1 1 ! .11 Ii "1 J " 

. .'n MMM, i i 

,'-!! iMMil i i M MM I 1 1 
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ill THEN 696 



MMi M 11 i 

MM i M j M I i i 

1,,'ti i ! " ..... 'MM 

: :, :i i i ,M 1 1 1 i ,' H i 

i. M i iM ,| ' "Mill i I 1 1 1 l |.'n ,' i i I i HK" ii 

''Mil MII'IM III 

'MM i i ||| " | |i i " MUM 'I'M 

" M IM'-I 1 "i I III II" 3 

"in i Mill i I 'i 

1 i i i i "I | M | | j i i i 1 i 

''in I I- MM "I Mil II "I '"llu 1 . Ml I H !! ERRSEIJ 1 ' BY 

1 n , i i i ,i i | I I i !' H I 

",'n ii M|' " ,| HIM 1 ! Ml ml I I- 1 ! 1 ; 11 5 

"mi MMM III 

MI M ||l " I |ii " | Ml M i M J I 

I MMM HIM 1 "I II MM MM I M M! I II US I BE>768> " 

In! M Mil'!!! | M 
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VfiLUE 1 



1839 D ISP "STORE VflLUES"? 

1640 INPUT Pi I 

1 050 IF fl $=--" NO" THEN 1 1 10 

1 6 D I S ! ::: ' ' ' F I !.... E- tt ( N "I" E F I !.... E M U S T B E > 5 1 2 ) ' ' 5 

11370 INPUT F8 

1080 IF F8>F9 THEN 1 100 

1090 REN 1 HI) 

I 100 STORE DRTFl F8? V 

II lO REN I HI) 

1120 PRINT "ft***;-*** MINIMUM I'M I'll ilMMi V 1 . i , 1; 
1 130 REN INITIFlLIZE SEQUEi li I' I ul -I r *': ! ' i 
1140 FOR I<=1 TO 20 
1 150 N[ !< >0 
1160 TC K > IE -I- 6 3 
1170 NEXT K 

1180 REN IHI.TIfiL.IZE SUM ITiBLE,, 
1190 FOR IO1 TO 60 
1200 JC lOO 
1210 SE K ]:==lE-i"S3 
1220 NEXT K 

12 30 REN IDENTIFY NODE OF ni-'l i.l II - ' 
1240 I) ISP "INITIFIL NODE IN I 'Ml 1 1"? 
1250 INPUT HO 
' 1260 N5-N0 
1270 JC HO >8[ HO :i'~0 

1 2 8 R E N 8 E R F? C H L I N K T ! :; i 1:3 L E F R !1 F R M N D E 
1290 K1--S0 

IS 00 Kl^Kl-H 

13 10 I l~INTa..ll<l !3 il/lOO) 
1320 IF 11 -HO THEN 13 50 
IS SO IF KK.1 THEN 1440 

1 3 4 GI ill T 1 3 ii:) 8 

1 3 5 O R E M i ::: i C C U N U L Fl T E T I M E F' R N R I G I N "i" :i T u ;1 r i D i: '. ^ * >- * * 3 

1360 TI^LI: KI ii s ::i 11*100 

1370 si=-si: N0 TWI; KI ::i 

1 3 8 R E M C H P i ::: l R E T E N P 8 U M 8 1 T T R E E T F! B I.... E S I m F R i H E :i "i" N D E 

1390 IF SI>SI::TI ::i THEN 1410 

1 4 @ GI 8 U B 1 6 1 

1410 IF KK.1 THEN 1440 

1420 REM CONTINUE LOOKING I nk " I I'liM" IHU'il . 

1 '4 3 61 G T 1 3 61 i!::i 

1440 REM CHECK SEQUENCE TFlHI I I i H ' II I 1 1 ' ! I 1 , , 

1450 FOR !<>! TO 26 

1460 IF N[|<]>0 THEN 1490 

1470 NEXT K 

1 4 S Gi 1" 1 5 2 @ 

1 4 9 R E M P R E P Fl R E T I D E N T I F Y N E X T :l F R U M :l N H E a 
H500 G08UB 1790 

1 5 1 GI i'J T 1 2 8 @ 

1520 REM GOTO HIM PFlTH PRINT OUT, 

1530 GOSUB 1940 

1 5 4 D I 8 P ' ' T H E R P Fl T H 8 F R N Fl N T H E R N D E ' ' !! 
1550 INPUT Fl$ 



1388 INPUT Fl* 

1590 IF fi*="YE8" THEN 138 

1 668 GOTO 2218 

1 6 1 i /i p F M * * * S U B 1 * * M Fi N 1 P U L Fl T 1 M ! :: T R E E T ! ::: i B !.,. E , 

1630 IF r SCT1 3== IE + 6 3 THEN 1700 

1630 M=l 

1640 FOR IO1 TO 20 

1658 IF SCT1 1-TCK3 THEN 1670 

1660 NEXT l< 

1670 SCT1 ]= SI 

1630 JCT1 3=11 

1630 CIO 'TO 1760 

1700 M ::::: 

1710 SCT1 3= SI 

1720 JCT1 3=11 

1730 FOR KXl TO 20 

1740 IF HE K 3=0 THEN 17 60 

1750 NEXT K 

1760 "IT. K 3 = SI 

1770 Nr.l<3 = Tl 

1730 RETURN 

1 7 9 pi F? F M * * * S U B 2 * * * I D E H T I F 1 Fl T 1 N F H E X T :1 F R N :! N U E 

1 8 R E M ! :: 1 N D M I N I M U M V Fl L U E I N S E Q U E N C E T Fi B !.... E , 

1310 T0=TC 13 

1320 M0 = M[ 1 3 
j :: ;: j/i [ ;;:: ;: } 

13413 FOR 'K=l TO 20 

1350 IF T0<T[i<] THEN 1330 

I860 T0:=T[I<3 

1370 M0='MCK3 

1330 I2*K 

1390 REN PRINT 1 1lLi/'-iS K5 Hf. K ::i5 TC K 3 

1300 NEXT K 

1910 T[ I23=lE+63 

1320 N[ 123=0 

1330 RETURN 

1 9 4 i/i R E N * ;>; * S I J B 3 * * # D E "I" E R N I H Fl T 1 N Fi N D F' R I H T I H G !' : M I i -S I M U M P R T H S 

1 9 5 H 1 S F I" E R M 1 N fl L N D E < 9 9 9 :::;: Fi L L H D E S ) ' ' 5 

1360 INPUT l<4 

1970 FOR ION1 TO 1 STEP 1 

1930 IF KttK.4 RUB K4tt999 THE hi 150 

1990 K2=K 

2000 Kl=0 

2010 K1=K1+1 

2020 HtKl 3=JLK2 3 

2030 IF J[K23=0 THEN 3060 

2040 K2=JCK23 

2050 GOTO 2010 



I 'I-.: IN ! I Mi: i mm MMlNIi! i'lli 1 1 ! i-'i i| 
!! ! MINI ,' I Ml 

! I II- 1 I : I I I Ml I '.Ml 1 i 
II I -IM I I I ' MINI ,'i ,'n 
I'!-' I II I "Mil liniH Mi ilniH i'll'i II 
i.'iMi ,'IJU 
I 'I- 1 MM I FIB 2 5 8 H[ l<3 ] 
Ml ',M ! .: 
I 'I- 1 MM I FIB 2 5 3 K 

iii ;M i 

il i -I '199 THEN 22S0 

HIM 1 "IINOTHER TIER MI NHL NODE 11 5 

1 1 II 'Ml n* 

M III "VIES 11 THEN 1956 

I'l MH-'ll 

i III! 
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The mission of the PACIFIC NORTHWEST FOREST 
AND RANGE EXPERIMENT STATION is to provide the 
knowledge, technology, and alternatives for present and 
future protection, management, and use of forest, range, and 
related environments. 

Within this overall mission, the Station conducts and 
stimulates research to facilitate and to accelerate progress 
toward the following goa!s: 

1. Providing safe and efficient technology for inventory, 
protection, and use of resources. 

2. Developing and evaluating alternative methods and 
levels of resource management. 

3. Achieving optimum sustained resource productivity 
consistent with maintaining a high quality forest 
environment. 

The area of research encompasses Oregon, Washington, 
Alaska, and, in some cases, California, Hawaii, the Western 
States, and the Nation. Results of the research are made 
available promptly. Project headquarters are at: 

Fairbanks, Alaska Portland, Oregon 

Juneau, Alaska Olympia, Washington 

Bend, Oregon Seattle, Washington 

Corvailis, Oregon Wenatchee, Washington 
La Grande, Oregon 



Mailing address: Pacific Northwest Forest and Range 

Experiment Station 
P.O. Box 3141 
Portland, Oregon 97208 



The FOREST SERVICE of the U.S. Department of Agriculture is dedicated 
to the principle of multiple u;-;e management of the Nation's forest resources 
for sustained yields of 'wood, _"wiier, forage, wildlife, and recreation. 
Through forestry research/ 'co-operation with' the States and private forest 
owners, and management of the Notiona! Forests and- National Grasslands, it 
strives as directed. by Gongr^s:; to provide mciesshgsy greater service to 
a growing Nation. '. . ~ ' ' ...v 

The U.S. Department -of Agriculture is an' Equal Opportunity Employer. 
Applicants for a!! Department programs wi!i. be yiven equal consideration 
without regard to race, color, sex or nation^! origin. 
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